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Abstract 

Two mechanisms of decoherence in ion traps are studied, specially related to the 
experiment [Kielpinski et ah, Science 291 (2001) 1013]. Statistical hypothesis are 
made about unknown variables and the expected behaviour of the visibility of the 
best experimental pattern is calculated for each mechanism. Data from the exper- 
iment are analyzed and show to be insufficient to distinguish between them. We 
suggest improvements which can do this with slight modifications in the present 
£h ' facilities. 
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I 

Ions traps are physical systems over which experimentalists have a large 
amount of precise control of the relevant degrees of freedom. This naturally 
makes these systems good candidates for exploring quantum phenomena and 
even for small scale quantum computation. Decoherence is the major barrier to 
be overtaken. One way of understanding decoherence is as a kind of quantum 
noise[l}. Its consequence is the transformation of pure states, which can be de- 
scribed by a state vector \ip), into mixed states to which only density operator 
description applies. Naturally, there is a huge set of decoherence sources, from 
spurious electromagnetic fields originated in the lab, passing through small os- 
cillations in parameter controls, up to background radiation and gravitational 
waves or whatever other kind of perturbation which forbids the system to be 
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completely isolated [2]. Usually they are all considered as an environment to 
which the system is coupled. 

In this work we give special attention to a paper by the Nist group [3], in which 
they showed how to produce a decoherence "free" qubit from a pair of ions 
trapped together. The idea of decoherence free subspaces can be classically 
understood, despite the fact of being a quantum phenomenon: if two classi- 
cal moving particles are subjected to exactly the same noise, their relative 
motion is unaffected. In fact, if one can determine the major source of deco- 
herence, and find two state vectors which are affected in the same way by free 
Hamiltonian and this decoherence source, any quantum information encoded 
in the subspace generated by these two vectors will be preserved from this 
decoherence source. As a cat and rat dispute, the other sources of decoherence 
will prevent this subspace to be completely decoherence free, but decoherence 
times will be substantially raised. That is what was obtained in two different 
situations in the cited experiment. 

We will start our analysis by considering that the major source of decoherence 
are oscillations in the magnetic fields used to split hyperfine structure of 9 Be + 
ions. More generally speaking, in situations like this, one can just consider the 
system to be subjected to a Hamiltonian H (£), which depends on a stochastic 
parameter £. Stochastic Hamiltonians generate stochastic evolutions, so pure 
states evolve to mixed ensembles in a quantum state diffusion picture [4]. If we 
consider a situation in which only the magnitude of the field is subjected to 
oscillations, then the eigenvectors of H (£) do not depend on £, but their eigen- 
values do. Let us consider a model like this for the ions: we are interested only 
in two internal states of each ion, so we consider the stochastic Hamiltonian 
to be {h = 1) 

= — 2 — zl H — 2 — ' ^' 



which can be rewritten as 

H (0 = ^# Ki + o z2 } + ^ {a zl - a z2 } , (2) 



where u m = (u\ + 002) /2 and Ud = (u\ — 002) /2. The form (2) makes evident 
the structure of the Hamiltonians H (£), including their dependence on £. In 
the ideal case, u\ =002, we get a doubly degenerate level with energy E — 0, 
and two other levels with energies ±u m . If the perturbations keep cu 1 = uj 2 we 
have a bidimensional decoherence free subspace generated by {|Tj) , lit)}- So 
the situation close to DFS is to have \u>d\ <C \w m \, and to consider initial states 
in this subspace. In this subspace, the term proportional to u m is ineffective, 
and decoherence is only originated by the variations in uj d {i.e.: only relative 
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variations are important, because they affect Born frequencies). 



Now we address the question: how does a state in a decoherence "free" sub- 
space decohere? Let us consider an initially pure state 



|lK0)) = a|U)+/W>- 



(3) 



For a fixed value £, the Hamiltonian (2) implies that at time t one has 

(t)) = ae-™^ | U> + Pe™*®* ||T) ■ 



(4) 



It is interesting to stress that in this model the vector state (4) always belongs 
to the decoherence "free" subspace, but due to its dependence on £, deco- 
herence emerges. Writing the state (4) as a density operator and taking the 
ensemble average over the stochastic parameter £, with weights p (£), one gets 



P(t) = 



o|2 



ka(3 



(5) 



where we defined the parameter 



k 
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(6) 



Let us stop for a moment to discuss states like (5). A Bloch vector picture 
is simple, specially in cylindrical coordinates: \a\ 2 — \f3\ 2 determines z coordi- 
nate, and 2 \ka(3*\ is the cylindrical radius (i.e.: the distance to the z axis). 
The most interesting case is \a\ 2 = \[3\ 2 = | (which corresponds to the exper- 
imental case), when \k\ directly gives the norm of Bloch vector, which can be 
interpreted as a direct measure of the purity of the state. In practice, the norm 
of a Bloch vector can be associated to the visibility of the best interferom- 
eter prepared with the state p. In interferometric experiments one naturally 
searches for the largest visibility [5], therefore we shall focus our attention on 
the factor k. 



Now let us discuss the factor k of the model here presented, eq. (6). As the only 
influence of £ we are considering is on the frequency cu^ (£) , we can consider as 
stochastic parameter a frequency v itself, with an ensemble weight p[y). So 
we have 

oo 

k= J e wt p(u)du, (7) 

— oo 
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where one can recognize the structure of Fourier transform of the stochastic 
weight p(v), so common in optics[6]. As k determines the visibility of the 
best interferometer prepared with state p, the data of ref. [3] must fit this 
expression for \k\. An interesting improvement on the experiment would be 
to independently monitor the magnetic field, and try to compare the fringe 
visibility and the Fourier transform of the statistical variations of the field. 
This would decide whether fluctuations of the magnetic field are the major 
decoherence source or not. As these fluctuations are not available, we can only 
speculate about them. One natural hypothesis is of Gaussian distribution. This 
would imply a Gaussian decay of \k\ with time t. Other distributions can also 
appear, but we have no a priori reasons for treating them. One last example 
is to suppose exponential decay for the visibility, what in this model would be 
consistent with Lorentzian stochastic weight 

r 

pMoc- -J— — . (8) 

[y - u ) + T 2 



In the cited experiment [3] , one state in the decoherence "free" subspace (DFS 
state) and one out of it (test state) were submitted to two distinct situations 
each: autonomous evolution, just subjected to natural noise, and an engineered 
reservoir, consisting of a dissonant laser with random intensity. The first sit- 
uation is supposed to be well modeled by the above discussion on variable 
magnetic field, but the second one deserves special attention. Experimental 
values allow one to consider only decoherence sources related to the presence 
of the dissonant laser. The most simple way to do it would be to consider 
the dispersive approximation to Jaynes-Cummings model for the interaction 
among one field mode and two two-level ions: 

HrfJC = H G + H int , (9) 



where 



H = |j 2 + cu /a ta, (10) 



where J = <j\ + a 2 , and 



Hi„i = 5f(j+a + J-a f ) 

wfl {aat (|Ti> (Til + IT2) (T2I) - ata (U\ + || 2 ) (| 2 |)} , (11) 

where f2 = g is a dipolar coupling constant, 5 = ^(uj — ujf) is the de- 
tuning, and dispersive approximation was applied, under the condition 5 2 ^> 



4 



g 2 (n + 1), n the number of photons. However, this can not mimic the exper- 
imental situation, unless we include somehow the randomization of the laser 
field. This can be done by considering an environment coupled to this field. As 
usual, let us consider a huge set of harmonic oscillators and bilinear coupling 
as a model to the environment. In this sense, the full Hamiltonian we must 
consider is 



The laser-environment part of this Hamiltonian can be rewritten in terms 
of normal modes, and this formally shows the equivalence between two ions 
coupled to a laser coupled to an environment and the same two ions directly 
coupled to a reservoir. To go further one needs to make some hypothesis on 
the coupling of the laser to the environment. White noise is a natural choice, 
and a Wigner-Weisskopf like approximation [7] implies the form of the coupling 
between ions and environmental normal modes. This is once again consistent 
with exponential decay of fringe contrast in interferometry. 

So, we considered two distinct sources of decoherence for a system like the one 
worked out in ref. [3]: one given by stochastic variations of the energy levels 
(probably due to variations of the applied magnetic field), and other given 
by the direct coupling to a "classical" field mode ("classical" in the sense of 
being strongly coupled to an environment). For both, we made statistically 
natural hypothesis, namely Gaussian fluctuations and white noise. The first 
model naturally gives Gaussian decay for the fringe contrast while the sec- 
ond implies exponential decay. We therefore strongly believe that when the 
laser is off, "natural" decoherence will fit a Gaussian, while with the laser 
on, "engineered" decoherence will be exponential. We take the data from the 
experiment, and used least square methods to obtain the best Gaussian and 
the best exponential to each set of data, and also used a kind of sieve to try 
to tell us if the points fit better a Gaussian or an exponential. We now discuss 
in detail the methods we used and the results we got. 

In fact, the function which we statistically treated was F (t) = InV (t), where 
V is the visibility of the fringe pattern. For each set of data we search for the 
least square curves of the form 



for which is important to say that both have the same number of parameters, 
what make the comparison fair. For each set of data, and each fitting family, 
we accumulate the vertical distances between the measured points and the 
best fitting curve and compare: for each set of data, the one which accumulate 



H = H dJC + ]T w,bfb, + Yl 9i (a f bj + bfa) . 



(12) 



F = at + b, 
F = At 2 + B, 



(13) 
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less distance is considered the best curve. For each initial state and laser 
situation the figure shows the experimental points, the best Gaussian (in fact, 
the best parabola for the logarithm) and the best exponential (resp. linear 
function). The accumulated distance is given in captions. The worked data 
are not enough to corroborate or deny our pre-conclusions about the form the 
coherence decays in each case, although a few more points could do the job. 
Also, as pointed out before, independently recording the field fluctuations and 
comparing its line shape to the visibility (with the laser off) could test the 
usually accepted viewpoint that the most important decoherence source are 
these field fluctuations. We expect that this work can stimulate more detailed 
experiments towards a better understanding of the mechanisms of decoherence 
on ion traps. 
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Figure captions 



Figure 1: For autonomous decoherence (laser off), the points and error bars are 
experimental [3], the curves are least square fits: (a) DFS state, y = —.803 — 
.00224t, with accumulated square distance (asd) .0095, and y = —.997 — .393 x 
10~ 5 t 2 , with asd .062; (b) Test state, y = -.109 - .00883t, with asd .037, and 
y = -.394 - .391 x 10~ 4 t 2 , with asd .0040. 

Figure 2: The same for engineered decoherence (laser on): (a) DFS state, 
y = -.874 - .00330t, with asd .0084, and y = -.884 - .159 x 10~ 3 t 2 , with asd 
.0083; (b) Test state, y = -.174-. 175t, with asd .084, andy = -.0110-.581t 2 , 
with asd .16. 
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